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Abstract 

We investigate the effects of spin-orbit interaction (SOI) and valley mixing on the transport and dynamical properties of 
a carbon nanotube (CNT) quantum dot in the Kondo regime. As these perturbations break the pseudo-spin symmetry 
in the CNT spectrum but preserve time-reversal symmetry, they induce a finite splitting A between formerly degenerate 
Kramers pairs. Correspondingly, a crossover from the SU(4) to the SU(2)-Kondo effect occurs as the strength of these 
symmetry breaking parameters is varied. Clear signatures of the crossover are discussed both at the level of the spectral 
function as well as of the conductance. In particular, we demonstrate numerically and support with scaling arguments, 
that the Kondo temperature scales inversely with the splitting A in the crossover regime. In presence of a finite magnetic 
field, time reversal symmetry is also broken. We investigate the effects of both parallel and perpendicular fields (with 
respect to the tube’s axis), and discuss the conditions under which Kondo revivals may be achieved. 
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1. Introduction 

The Kondo effect [1] is a hallmark of strongly correlated 
electron physics. Its observation in quantum dot set-ups is 
ubiquitous and reveals precious information on the under¬ 
lying symmetries of the quantum dot system and on the 
corresponding degeneracies of its spectrum. Specifically, 
electrons in carbon nanotubes (CNTs) possess a spin and 
a pseudo-spin degree of freedom [2], the latter originat¬ 
ing from the presence of two inequivalent Dirac points in 
the underlying graphene hexagonal lattice. In the absence 
of spin-orbit interaction, and considering only transverse 
quantization, the CNT’s Hamiltonian is invariant under 
time-reversal and pseudo-spin reversal symmetries, and 
thus a quadruplet of degenerate levels is associated with a 
given longitudinal momentum. In this case, the four-fold 
degeneracy may lead to the occurrence of the so called 
SU(4)-Kondo effect at low temperatures [3-7]. In order to 
see this exotic Kondo resonance it is important, however, 
that both spin and pseudo-spin quantum numbers be con¬ 
served during tunneling (or reflection), as a mixing of these 
degrees of freedom can result in a more conventional SU(2) 
Kondo effect [8]. For CNTs devices where parts of the tube 
act as leads (see Fig. 1), such a situation can be realized, 
and the peculiar features associated to the presence of both 
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spin and orbital degrees of freedom can be probed in finite 
magnetic fields [6, 9]. Recently, SU(4) Kondo physics, also 
originating from coupled spin and orbital degrees of free¬ 
dom, could be engineered in double-quantum dot based 
devices [10]. 

In a more realistic description of a CNT though, pseudo¬ 
spin symmetry breaking contributions, like the curvature 
induced spin-orbit interaction (SOI) [11, 12] or valley mix¬ 
ing due to scattering off the boundaries [13, 14] or to disor¬ 
der [12], should be included. As a consequence, the four¬ 
fold degeneracy is broken and, for a given value of the 
longitudinal momentum, the spectrum of an isolated CNT 
quantum dot consists of two pairs of degenerate Kramers 
pairs, with splitting provided by the combined effects of 
SOI and valley mixing [15]. In this situation, upon in¬ 
creasing the SOI strength or the valley mixing, a crossover 
from the SU(4)-Kondo state involving both Kramers pairs, 
to the more standard SU(2) Kondo regime is expected [16- 
18]- 

Despite the considerable amount of experiments report¬ 
ing Kondo behavior in CNTs, [6, 9, 18-24], the combined 
effect of SOI, valley mixing and the impact of applied mag¬ 
netic fields on the SU(4) to SU(2) crossover have only been 
addressed within a field theoretical effective Keldysh ac¬ 
tion approach [18]. A numerically exact investigation of 
the highly intricate crossover is thus very desirable. 

In this work we study the dynamical and linear trans¬ 
port properties of CNT-based Kondo quantum dots by 
means of the Density Matrix-Numerical Renormalization 
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Group (DM-NRG) method [16, 25-27]. We focus on the 
SU(4) to SU(2) crossover induced by finite SOI and valley 
mixing, and study the influence of magnetic fields paral¬ 
lel or perpendicular to the CNT’s axis. At zero magnetic 
held, it is not possible to distinguish at the level of the 
spectral function or of the linear conductance among the 
two symmetry breaking effects. Here, what matters is the 
amplitude A of the total inter-Kramers splitting. We de¬ 
termine here the energy scales for the cross-over region and 
demonstrate that the Kondo temperature scales inversely 
with A, in agreement with previous analytical predictions 
[28, 29]. 

At finite fields, the behavior of the Kondo resonance is 
strongly influenced by the relative strength of the SOI and 
valley mixing contributions, as well as by the direction of 
the applied held. In fact, a major effect of the curvature in¬ 
duced SOI is to set as spin quantization axis the tube’s axis 
and to lock spin and valley degrees of freedom [11], Valley 
mixing instead does not act on the spin degree of free¬ 
dom, but induces a rotation in valley space [13, 15]. Thus, 
in parallel held the spin is still a good quantum number. 
Magnetic held induced Aharnov-Bohm contributions dom¬ 
inate over Zeeman effects at small helds, due to the large 
orbital moment of the nanotubes [ 2 ] , which enables one to 
clearly resolve the splitting [23] and rejoining [15, 18] of a 
Kramers pair also at low helds. In perpendicular helds, in 
contrast, the spin is no longer a good quantum number; 
rather the simultaneous presence of SOI and valley mixing 
implies a full entanglement of orbital and spin degrees of 
freedom. In this case it is convenient to classify virtual 
Kondo transitions in terms of the discrete operations re¬ 
lated to the time reversal and pseudo-spin reversal opera¬ 
tors [18]. These considerations are nicely confirmed by our 
simulations and rehected, in particular, at the level of the 
linear conductance in the occurrence of Kondo revivals: at 
specihc values of the magnetic held, which depend on the 
held direction, valley mixing strength and on the number 
of charges trapped in the dot, a Kondo resonance can be 
restored near avoided level crossings. 

The paper is organized as follows. We present our model 
Hamiltonian for CNTs in Sec. 2.1, and give a brief analy¬ 
sis of the symmetries of the system in Sec. 2.2. Dynamical 
and transport properties are discussed in Secs. 3 and 4. 
We extend our discussion in Sec. 5 by including the effects 
of an applied magnetic held parallel (Sec. 5.1) or perpen¬ 
dicular (Sec. 5.2) to the CNT’s axis. Our conclusions are 
summarized in Sec. 6 . 

2. Theoretical framework 

2.1. Model Hamiltonian 

The setup we consider consists of a CNT quantum dot 
coupled to two external leads (see the sketch in Fig. 1). 
We focus on a single longitudinal mode (also known as 
“shell”), and correspondingly, describe the CNT by an ex¬ 
tended Anderson impurity model [1, 15], consisting of a 



Figure 1: Sketch of the setup. The carbon nanotube is coupled to 
two external contacts, source and drain. The quantum dot formed 
in the nanotube is indicated by the shaded area. The electrostatic 
state of the dot is capacitively controlled by a back-gate. 

pair of interacting Kramers doublets. We denote by Ej 
the energies of the four levels (j = {1,2, 3,4}), and by 
hj = djdj their occupation. In what follows, we shall refer 
to this basis as the Kramers basis (see Fig. 2(a)). Each 
of the four levels can accommodate one electron and, with 
a good approximation, these electrons interact with each 
other through a strong and level-independent on-site in¬ 
teraction U. In this basis, the CNT Hamiltonian takes the 
form 

4 4 

-ffcNT = £ 3 + U ^2 r b r V ' (!) 

3 — 1 3<j' 

In the absence of the spin-orbit interaction and valley 
mixing, Ago = 0 and Akk’ = 0, the CNT’s Hamiltonian is 
invariant under time-reversal and valley-reversal [2] . These 
operations are represented by the two antiunitary opera¬ 
tors T and V, respectively [18], and yield a fourfold de¬ 
generate spectrum of the CNT, Ej = £<}. Correspondingly, 
the CNT Hamiltonian is SU(4) invariant. In what fol¬ 
lows, we shall label states such that (1,2) and (3,4) form 
Kramers pairs, while (1,4) and (2,3) are pairs associated 
with the V symmetry. Notice that a third unitary opera¬ 
tor C = V'T~ 1 linking the remaining pairs (1,3) and (2,4) 
can also be constructed from T and V (see Fig. 2(b)). A fi¬ 
nite A = \/Ag 0 + A( <k , breaks the V symmetry and, cor¬ 
respondingly, also the SU(4) symmetry (see Appendix A 
for details on how these states and the symmetry oper¬ 
ations are constructed). Since time-reversal symmetry is 
preserved, the on-site energies remain twofold degenerate, 
Ei = £2 = £d + A/2 and £3 = £4 = £d~ A/2 (see Fig. 2(a)). 
Notice that a finite A plays the same role as a magnetic 
held on the V - and C-pairs, such that conjugation rela¬ 
tions among energy levels exist: £i(A) = £ 4 (—A), and 
similarly for the other couples. 

In a CNT, with a good approximation, each nanotube 
level couples to independent channels in the leads, and 
their tunnel coupling can thus be described by the Hamil¬ 
tonian 

4un = ^ V 3 V~P0 j de + hx - ' ( 2 ) 
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Here, instead of the original left/right operators cA L ^ R for 
the leads, we introduced the symmetric and antisymmetric 
combinations 

A = / cos7j shayA /c e y L \ 

\k;i) \-siii7j COSJj) \C e j, R J ’ 

and the corresponding effective tunneling amplitude 
Vj = yjVj 2 l + Vj 2 r and asymmetry parameter 7 y = 

arctanV/ r/V/l- The leads are assumed to be non¬ 
interacting with a constant density of states per flavor 
p(ta) = po = 1/2 W, and a bandwidth 2 W. They are 
described by the Hamiltonian 1 

w 

ITlead = 'y ' J de £ {^ e jd £ j + bljbej'j . (3) 
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Notice that only the a £ j channel couples to the dot, while 
channel b £ j remains completely decoupled in equilibrium. 
The total Hamiltonian 


H = HcNT + H tun + Hlead (4) 

captures the essential physics of our set-up and, under 
equilibrium conditions, can be solved using Wilson’s NRG 
method [25]. 


2.2. Global symmetries 

Let us now discuss the continuous symmetries of the 
Hamiltonian (4). These symmetries are extremely useful, 
since they allow for an efficient numerical treatment of the 
problem. Throughout this paper, we shall focus on the 
simplest but physically relevant case of 


Vj =V 7j — 7 • 


In this case, for A = 0, the total SU(4)-spin operator 


jSU(4) 






(5) 


commutes with the Hamiltonian (4), and the SU(4) sym¬ 
metrical Anderson model [7] is recovered. The A’s above 
denote the 15 generalized Gell Mann matrices or some 
other set of matrices defining the SU(4) representation, 
and the operators (5) satisfy the SU(4) Lie algebra. 

Finite inter-valley scattering or spin-orbit field imply 
A ^ 0, and break the SU(4) symmetry down to SU(2)(g> 
SU(2). The latter are generated by the usual SU(2) spin 
operators 


Me 
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Quasiparticle operators are normalized to satisfy {a £ j, a £ /j/} = 



Figure 2: a) Energy level scheme involving two Kramers doublets 
separated by the energy A. (b) Connection among the energy levels 
established by the symmetry operations T, V and C. 


acting on the two Kramers doublets n = (1,2) and k = 
(3,4) (see Fig. 2). Here, <x = ( a x , a yi a z ) is the regular vec¬ 
tor of the Pauli matrices. In addition to these SU(2) sym¬ 
metries, the total charge is also conserved in each Kramers 
’channel’ 

= \ E ( d 1 ^ -1/2+ f de : aVd ej ■) , (7) 

with : ... : referring to normal ordering. 


3. Transport properties 


In this section we present the results for the spectral 
functions Aj (oj) of the operators d)- and evaluate the con¬ 
ductance across the dot under equilibrium conditions. The 
linear conductance can be computed directly within the 
NRG and is related to the equilibrium spectral function of 
the operators d/-. It reads 


G{T) = 



— OO 


df(u,T) \ 
duj J 


Aj (w), 


( 8 ) 


with F ; = nV^po the usual broadening parameter, and 
aj = 4 tan 7 j /(1 + tanyj ) 2 the asymmetry prefactor, 
which depends on the source and drain couplings and 
is in general smaller than one. In Eq. ( 8 ), f(ui,T) = 
(l+exp(w/T )) -1 (unit fcg = 1) is the Fermi-Dirac distribu¬ 
tion function and Aj (uj) denotes the equilibrium spectral 
function of the j-th dot level, 

A j M = - ( w )}. ( 9 ) 

with Gfj, (u>) being the Fourier transform of the retarded 
Green’s function GV, (t) = —iQ(t)({dj(t ) ,dj,(0)}). To 
compute Aj (w), we used the open-access Budapest DM- 
NRG code [30], that explicitly uses the symmetries of the 
system. As discussed above, for A = 0 the system exhibits 
SU(4) symmetry. Therefore, by tuning the parameter A 
from 0 to some large value A >> W, for a singly occupied 
longitudinal level, we can follow the crossover from the the 
SU(4)-Kondo fixed point to the SU(2)®SU(2) one. In the 
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Figure 3: The total spectral function Atot (tu) as a function of fre¬ 
quency for different values of A (color scale in units of 4 = 

Tk (A = 0) ~ 0.000429r). Other parameters were fixed to U = W, 

= —U/2 (iV g = 1), r = U/50 and T = 0. Inset: The evolution of 
the Kondo peak for different values of A. 

next two subsections we shall study the manifestation of 
the crossover at the level of the spectral functions and of 
the linear conductance, respectively. 

3.1. Spectral functions 

The impact of the splitting A on the Kondo reso¬ 
nance is demonstrated in Figs. 3 and 4, where we dis¬ 
play the total spectral function, A to t(oj) = JTA,-(w), 
for two specific values of £d and a relatively large ra¬ 
tio U/Y = 50. For A = 0 the occupation of the CNT 
levels is controlled by the ’dimensionless gate voltage’, 
N g = (—e d + C//2 + A) / ([/ + A/2), taking on integer val¬ 
ues N g = k = 1,2,3 just in the middle of the Coulomb 
blockade valleys with k particles on the CNT. 

Fig. 3 displays the crossover from the SU(4) to the SU(2) 
regime for the case N g = 1 and several values of A. By 
particle-hole symmetry, the spectral functions for N g = 3 
are the mirror images of the N g = 1 spectral functions, 
and we do not discuss them in detail. In the limit A = 0, 
an SU(4) Kondo resonance arises due to quantum fluc¬ 
tuations of the ground state quadruplet. As expected [1], 
this resonance is pinned asymmetrically to the Fermi level, 
uj = 0. As soon as the splitting A becomes comparable to 
(extracted form the half width at half maximum of 
Atot(w)), the spectral function maximum lowers and tends 
to be symmetrical around the Fermi energy. At the same 
time, two satellite peaks emerge at approximately ±A. 
These satellite peaks correspond to “electron-hole” exci¬ 
tations between the two Kramers pairs depicted in Fig. 2. 
Notice that the value of the spectral function does not 
change at the Fermi energy as the SU(4) resonance grad¬ 
ually turns into a symmetrical SU(2) Kondo resonance. 
This implies that for N g = 1 the T = 0 temperature con¬ 
ductance is not suppressed by breaking the SU(4) symme¬ 
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Figure 4: The T = 0 temperature total spectral function Atot (^) as 
a function of frequency for different values of A and at = —3U/2 
(N g = 2). We have used U = W, and F = U/ 50. (Color scale in 
units of — 7 ^ k (A = 0) ~ 0.00031ir.) Inset: Evolution of the 

Kondo peak on a logarithmic scale for different values of A. 

try. However, as shown in the inset, the Kondo tempera- 
ture is strongly reduced for A 3> T^}' <4> . 

The situation is dramatically different for N g = 2 
(£d = —3(7/2). At this value of £d there are two electrons 
on the CNT longitudinal shell. The Hamiltonian exhibits 
electron-hole symmetry, and the total spectral function 
is symmetrical for any value of A. As a consequence of 
Friedel sum rule [31], the value of A tot (w = 0) is twice 
as large as it was for N g = 1. For N g = 2, however, a 
splitting A Tf (4) eliminates the ground state degener¬ 
acy of the isolated CNT, completely suppresses the Kondo 
resonance, and leads to the emergence of a ’pseudogap’ of 
width ~ A in A to t(w). 

3.2. Linear conductance 

The crossover features in the spectral function are also 
reflected in transport characteristics. In Fig. 5, we present 
the T = 0 temperature linear conductance as a function 
of N g for several values of the splitting A. Similar to 
the spectral function, as a manifestation of the two elec¬ 
tron SU(4)-Kondo state, the linear conductance also ac¬ 
quires its maximal value at the particle-hole symmetric 
point N g = 2 in the SU(4) symmetrical case, A = 0 [17]. 
This large conductance is, however, sensitive to A and is 
quickly suppressed for A as a consequence of 

the pseudogap appearing in the spectral function. 

In contrast, for N g = 1 and N g = 3 a large A does 
not destroy the conductance at T = 0 temperature, which 
remains close to 2 e 2 /h. Notice, however, that the origin 
of the conductance is different in the limits A = 0 and 
A X^ U( ' 4 \ In the SU(4) limit, A —> 0, incident conduc¬ 
tion electrons pass through the quantum dot with proba¬ 
bility 1/2 in all four channels. In contrast, for A —> oo the 
two empty levels and the corresponding channels do not 
conduct at all, while the other two have a perfect, unitary 
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Figure 5: Linear conductance as a function of Ng for different values 
of the splitting A (T^ U ^ = Tk (A = 0, Ng = 1) = 0.000429r). The 
parameters are U = W, T = U /50 and T = 0. 


conductance due to the residual SU(2) Kondo effect. Since 
in these channels electrons pass through the dot with prob¬ 
ability one, in this A —> oo limit the linear conductance 
becomes noiseless [1]. 

Another useful way to visualize the crossover between 
the SU(4) and SU(2) regimes is to consider the tempera¬ 
ture dependence of the linear conductance, shown for N g = 
1 in Fig. 6. To explore universal scaling, we have rescaled 
the temperature by the Kondo temperature T^ CM (A), de¬ 
fined as the temperature at which the conductance is re¬ 
duced to half of its T = 0 temperature value [1]. For A 
much smaller than Tj!? CM (A), conductance curves show 
SU(4) universality and lie on the top of each other. As 
soon as A becomes comparable with T^ CM (A), however, 
universality is lost, the Kondo temperature lowers and a 
peak emerges at approximately A. For very large values 
of A/Ijb CM (A), the curves become universal again, but 
they now follow a simple SU(2)-scaling. The two univer¬ 
sal curves for A = 0 and A > U are associated with the 
SU(4) and SU(2) fixed points, and can be derived from the 



T/T k ( A) 


corresponding Kondo models (see Appendix B for further 
details). 


4. Kondo temperature 


The Kondo temperature is the most basic energy scale 
that characterizes the correlated Kondo state. As re¬ 
marked already, in the regions N g « 1 and N g « 3 it is 
dramatically suppressed for large Kramers pair splittings 
A > 7| U(4) . This suppressed SU(2) Kondo temperature, 
can be estimated by carrying out a two-step RG 
procedure [32]. 

The Kondo temperature has first been defined in the 
context of a simple spin exchange Hamiltonian of the form 
Hint = J S • s, describing the interaction of a magnetic 
moment S at the origin with the local spin density s of 
the conduction electrons, with J > 0 an antiferromagnetic 
exchange coupling. Later on, it was realized that the sym¬ 
metry of the exchange Hamiltonian plays a crucial role 
and affects the Kondo temperature. The SU(N)-Kondo 
model [33], in particular, 2 yields a Kondo temperature 
X^U(N) ~ \Y e -2/Np 0 J^ implying that a larger N apprecia¬ 
bly enhances Tk. This behavior can be obtained by carry¬ 
ing out a renormalization group analysis for the effective 
exchange coupling (vertex function) describing scattering 
processes at energy e . To leading logarithmic order, the 
dimensionless effective exchange coupling, j(e) = poJ e g, 
satisfies the RG equation 


d l = 1 l ,2 

dx 2 ' ’ 


( 10 ) 


with x = In {W/e) the scaling parameter, and diverges log¬ 
arithmically at the Kondo scale e = T^ U( - N \ 3 

Let us now focus on the Kondo regime (T <C U) of a 
CNT quantum dot with N g = 1 electrons trapped inside a 
shell. In this regime, assuming further A <C U, the interac¬ 
tion of the electrons with the nanotube can be described by 
an exchange Hamiltonian of almost perfect SU(4) symme¬ 
try [7]. The structure of the vertex function (effective ex¬ 
change amplitude), however, depends on the energy of the 
electrons scattered. Electrons of very high energy, e> A, 
can induce transitions between all four levels, and expe¬ 
rience an SU(4) exchange process. Correspondingly their 
scattering amplitude satisfies (10) with N = 4. Electrons 
of energy e -C A, however, can only flip the states within 
the lower Kramers doublet, and their scattering amplitude 
obeys (10) with N = 2. Matching the N = 2 and N = 4 
solutions of Eq. (10) at x = \n(W/ A) we thus obtain 


j(e < A) 


In 


A/(T] 


K 


Figure 6: DM-NRG result for the linear conductance as a function 
of T/Tk for different values of the splitting A. The parameters are 
U = W , r = U/ 50, T = 0 and = —U/2. Here Tk(A) stands for 
TJCM(A) defined in the main text. 


2 The SU(N) Hamiltonian is defined in terms of an iV-fold degen¬ 
erate level as H [ nt = (J/2) X , with X the exchange operator. 

3 Both e and W can be considered as scaling variables here. 
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The energy where the perturbative RG breaks down can be 
identified as the Kondo temperature Xk, and is expressed 
as 

rK(A)a^e-^, (11) 

predicting a ~ 1/A decay of the Kondo temperature for 
A < U. Notice that in this expression, we replaced the 
bandwidth by U, which serves as an upper cutoff for spin 
exchange processes within the Anderson model. 



Figure 7: Evolution of the Kondo temperature as a function of A. 
The Kondo temperature was extracted from the NRG data in two 
different ways: i ) as the half width at half maximum of the total 
spectral function (HWHM), ii) half conductance maximum (HCM). 
The other parameters were fixed to: U = W, z,\ = —U/2 [Ng = 1), 
r = U/50 and T = 0. 

As we now demonstrate, the analytical expression (11) 
is in good agreement with our NRG calculations. The 
definition of Xk is not unique, we have therefore extracted 
it from the NRG data in two different ways: in addition 
to the half conductance scale X^ CM (A) defined earlier, we 
also introduced the Kondo scale X/? WHM (A) as the half 
width at half maximum of the spectral function for the dj 
operators. Both Kondo temperatures are shown in Fig. 7, 
where three regimes can be delineated: For A < X^ U( ' 4 \ 
the system is governed by SU(4) physics, and Xk agrees 
with the SU(4) Kondo temperature. This is followed by 
a ’crossover region’, < A < U, where an SU(4)—» 

SU(2) crossover takes place as a function of temperature or 
energy, but below A an SU(2) Kondo state emerges with 
a suppressed Kondo temperature given by (11). Finally, 
for A > U only one Kramers doublet remains active, and 
an SU(2) Kondo behavior appears at all scales with a Xk 
independent of A. 

5. Transport in a finite magnetic field 

A finite magnetic field turns out to be crucial to dis¬ 
entangle the effects of SOI and valley mixing. As shown 
below, both the spectral function and the linear conduc¬ 
tance display qualitatively different features depending on 
the direction of the applied magnetic field. 

In finite magnetic field the N g = 1 (center of the 
first Coulomb valley) condition becomes = —U/2 — 


(£ 4 (B||) + £ 3 (-B||))/2, whereas N g = 2 (particle-hole sym¬ 
metric point) remains £d = — 317/2. 

5.1. Parallel magnetic field 

Spectrum. The energy levels of the CNT Hamilto¬ 
nian, Eq. (1), easily follow from the discussion in Ap¬ 
pendix A.2. They read 

1 1 /—“- 2 

El,4 = Ed + ~^ 9 sB\\ 2 V ^KK’ "k (Aso + 2g , or b-B||) , 

(12a) 

£2,3 = e d — i - y A|- k , + (Ago — 2g or b-B||) , 

(12b) 

where g s (<?orb) is the Lande spin(orbital) g-factor and B\\ 
the amplitude of the parallel magnetic field. 

The resulting spectrum for occupancies with N = 1,2,3 
of the CNT is shown in Fig. 8. It is interesting to notice 
the crossing of the states |3) and |4) in the sector N = 1 



-1.5 -1 -0.5 0 0.5 1 1.5 

BJB S 


Figure 8: Parallel magnetic field dependence of the energy levels 
for (a) N = 3, (b) TV = 2 and (c) TV = 1 Hilbert space sectors. 
Black circles indicate ground-state levels crossing. The parameters 
are U = W, F = 0.02VF, = —U/2, g s = 2, g OT b = 1.83<? s and 

Ark 1 = Aso/2 = 5r. They are chosen such that the inequality (13) 
is fulfilled. 
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Figure 9: Evolution of the single particle energy levels in paral¬ 
lel magnetic field. The parameters are T = 0.02VF, U = W, 
e d = - U/2, A S q = 5r, g s = 2, ^ orb = 1.83 g s and A K k’ = 
lOAso yj ^orb/^s — 1- F° r this choice of parameters the inequal¬ 
ity (13) is not satisfied and no level crossing is possible at finite 
magnetic field. Here, B e > 0 corresponds to the value of B n such 
that the state |2) shows a minimum. 


Focusing on the N = 2 sector of the Hilbert space (see 
Fig. 8(b)), one notices that again an avoided crossing oc¬ 
curs at B OI b = Aso/2fforb) which becomes an exact cross¬ 
ing for Aki<’ = 0. In this special case the two states |34) 
and 124) are degenerate. Finally, no ground state crossing 
is observed for N = 3 (see Fig. 8(a)). 

Using these considerations, we can build up the ground 
state configuration of the system as a function of the 
applied parallel magnetic held and of the orbital energy 
£d (see Fig. 10). By inspecting where a crossing occurs 
(dashed lines in Fig. 10), it is possible to capture the so- 
called “Kondo revivals” [17], where the Kondo effect is 
restored for some specific values of the applied magnetic 
held. As Fig. 10 reveals, revivals are expected around 
B\\ = B OI b for Akk’ — 0 in the N = 2 valley, and around 
B\\ = B s , if the inequality (13) is fulfilled, in the N = 1 
valley. 

These considerations are rehected in the behavior of the 
spectral function and of the linear conductance. 


(circles in Fig. 8(c)). It can occur if 


l A KK'| / 4 fforb 

lAsol V 3 S 2 

at a magnetic held value given by 


B s 



A 2 

a KK’ 


4 5orb - 9\ 


2 ‘ 


(13) 


(14) 


Thus, the ground states in the one-particle sector of the 
Hilbert space switch from |4) to |3) for B\\ > B s . Moreover, 
an avoided crossing occurs at 73 1 = B OI b = Aso/(2fforb)- 
If Eq. (13) is not satished, the crossing of the two states 
in the N = 1 sector cannot be achieved, as it is shown in 
Fig. 9. 


Spectral function. The spectral function is shown for sev¬ 
eral values of the parallel magnetic held in Fig. 11. In 
the absence of the magnetic held (black curve) it shows a 
Kondo peak at the Fermi level and two satellite peaks lo¬ 
cated at w ~ ±A = ±(e 1 ( 2 )(0) —£3(4)(0)). As we switch on 
the magnetic held, the Kondo peak lowers and then splits 
into two resonances located at w ~ ±(e 3 (13||) — e 4 (i3||)), 
due to processes representing quantum fluctuations within 
the lowest Kramers pair. Increasing the magnetic held, 
the peaks merge again into one at B» = B s (violet curve) 
and the spectral function recovers the unitary value at 
oj = 0 (see inset of Fig. 11). For larger values of B\\ the 
Kondo peak splits again. Regarding the satellite peaks, 
they lower, broaden and shift towards higher energies. We 
hnd that at B y = B s the leftmost one is due to a reso¬ 
nance in the third and fourth components (A 3 (u;) = A^(uj), 



B \\/ B S 



co/A 


Figure 10: Ground state configuration as a function of the energy 
and the applied parallel magnetic field. Grey dashed lines indicate 
degeneracies. The black dashed line indicates B n = B orb . Param¬ 
eters are, as in Fig. 8, U = W, T = 0.02 W, = —U/ 2, g s = 2, 

Sorb = 1.830s and A K k’ = A S o/2 = 5r. 


Figure 11: Close-up image of the spectral function as a function 
of the frequency for Akk’ = 0 and several values of the parallel 
magnetic field in the middle of the Coulomb valley with N g = 1 at 
T = 0. Inset: Full spectral function. The unitary limit is reached 
for B ||=0 while the full Kondo revival is obtained at = B s . 
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£3 = £4 for By = B s ) of the spectral function. On the other 
hand, the rightmost peak has an internal substructure be¬ 
ing the sum of a maximum in the first, and second, 

A 2 (a;), components. Increasing furthermore the magnetic 
field (By > B s ) the satellites peaks lower and broaden, the 
leftmost being determined essentially by the component 
corresponding to the ground state, A 3 (e). Despite the fact 
that the energy separation, e 2 (By) — £i(By), of the excited 
states |1) and |2) is monotonically increasing, the total 
spectral function does not show a split of the correspond¬ 
ing resonances in this range of investigated parallel mag¬ 
netic field values (0 < By < 2 B s with B s ~ 90Tk(A)). Our 
results on the impact of a parallel magnetic field on the 
spectral function of a CNT are consistent with the NRG 
analysis performed in Ref. [8] when Aso = Akk’ = 0. 




5||/5 0 rb 


Figure 12: Linear conductance as a function of the parallel magnetic 
field (a) in the middle of the Coulomb valley, JV g = 1, and (b) at 
the particle-hole symmetric point, N g = 2 at T = 0. Ago and Akk’ 
were varied such that A/T^ l '' ( ''' remained constant. 

Linear Conductance. In Fig. 12 we show the linear con¬ 
ductance of the system as a function of the applied parallel 
magnetic field and for a fixed ratio of A/T^ U ^ 4 ' 1 ~ 12. In 
Fig. 12(a) the conductance shows, as expected, the Kondo 
revival at By = B s . Moreover, the width of the Kondo 
peak is proportional to Tk(A). By increasing the valley 
mixing term Akk’ the resonance shifts towards By = 0 
and, in case Akk'/^so does not fulfill Eq. (13), the Kondo 
revival disappears (black solid line). 


In the N g = 2 valley, Fig. 12(b), the valley mixing term 
acts slightly differently. For Akk’ = 0 the SU(4)-Kondo 
effect in By = 0 is essentially suppressed for this set of 
parameters (A/T^ 1 ^ 4 -’ ~ 12) and hence a Kondo peak is 
present only at By = B or b- Switching on the valley mixing, 
since B or b is not Akk’ dependent, results in the suppres¬ 
sion of this resonance without any shift. We notice that 
Kondo revivals have indeed been seen in experiments on 
CNT-dots in parallel fields [15]. 

5.2. Perpendicular magnetic field 

Spectrum. A magnetic filed perpendicular to the CNT axis 
couples only to the spin degree of freedom. As it follows 
from Appendix A.3, the energy levels of the CNT Hamilto¬ 
nian (1) in the presence of a perpendicular magnetic field 
read 

£ 1,4 = £d ± 2 \I ^io + (Akk’ + 9sB ±) 2 ) (15a) 

£ 2.3 = £d ± ^\j Ag 0 + (A K k’ - g s B_ l) 2 , (15b) 

where B± is the amplitude of the perpendicular magnetic 
field. The corresponding spectrum is shown in Fig. 13. 
It shows clear qualitative differences with respect to the 



B J_/B kk’ 


Figure 13: Perpendicular magnetic field dependence of the energy 
levels for (a) N = 3, (b) N = 2 and (c) N = 1 particle number. 
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parallel case. For example, no ground state crossing is 
observed for N = 1 and N = 3: the distance between 
the lowest energy levels first increases and then saturates 
to the valley mixing strength Aki<’- An avoided cross¬ 
ing between excited states occurs for B± = ±Bkk’i where 
B KK ’ = Akk’/Ss- On the other hand, in the N = 2 sec¬ 
tor, the avoided crossing at B± = ii?KK’ is between the 
ground state 134) and the state 124). Hence an exact cross¬ 
ing occurs for Aso = 0. Comparing Figs. 13 and 8 we 
indeed see that the SOI and the valley mixing exchange 
their role changing from parallel to perpendicular mag¬ 
netic fields. Following the same argument of the previous 
section, this revels that the Kondo revival can be achieved 
only in the N s = 2 valley for B± = Bkk ’■ 



w/A 


Figure 14: Close-up structure of the total spectral function as a func¬ 
tion of the frequency for several values of the perpendicular magnetic 
field, for 7V g = 1 and T = 0. Inset: Full spectral function, showing 
that the unitary limit is only obtained for B± = 0. 


Spectral function. We analyze the behavior of the total 
spectral function at N g = 1 for several values of the per¬ 
pendicular magnetic field in Fig. 14. Since the Kondo re¬ 
vival is not possible in the first valley, the Kondo peak 
monotonically splits into two sub-peaks located at u> = 
±(£ 3 ( 1 ?j_) — £a(B±)). Furthermore, two satellite peaks are 
visible. The rightmost resonance is the sum of two peaks 
in Ai(e) and A 2 (e), whereas the leftmost one is essentially 
determined by the fourth component, A 4 (e); correspond¬ 
ing to the ground state level. As in the parallel case, no 
splitting of the satellites is observed in the investigated 
magnetic field range. For B± > Bkk’/2 the central Kondo 
peaks merge with the lateral satellites. 

Linear conductance. In Fig. 15 we show the linear conduc¬ 
tance as a function of the applied perpendicular magnetic 
field. A perfect Kondo revival is obtained only if Aso = 0; 
for Aso 7 ^ 0 the Kondo resonance is suppressed because 
the degeneracy between the sates |34) and |24) is lifted 
and the fluctuations of the quantum numbers referring to 
these states become energetically unfavorable. 



B ±/ B KK’ 


Figure 15: Linear conductance as a function of the perpendicular 
magnetic field at the particle-hole symmetric point, 7V g = 2 at T = 0. 

6. Conclusion 

In this work, we have investigated the transport and 
spectral properties of CNTs in the Kondo regime. To do 
that, we have constructed a model that, for the sake of 
simplicity, assumes spin and valley quantum number con¬ 
servation during tunneling processes, but accounts for the 
combined effects of spin-orbit interaction, valley mixing, 
and electron-electron interaction, U : and finite magnetic 
field. 

First, we carried out a detailed study of the cross-over 
between the SU(4) and SU(2)-I<ondo regime by analyzing 
the dependence of the spectral functions and the linear 
conductance on the Kramers pair splitting, A. We have 
shown by means of DM-NRG computations, in particular, 
that in the Kondo regime of a singly occupied CNT longi¬ 
tudinal mode a universal SU(4) conductance is displayed 
for A < while for A 7g> Tj® U(4 ^ the conductance 

and the spectral functions display characteristic SU(2) be¬ 
havior with a strongly reduced Kondo temperature. In 
the intermediate regime < A < 17 we observed 

a 1/A dependence of the Kondo temperature. This be¬ 
havior can be explained in terms of simple scaling argu¬ 
ments, and is also in agreement with exact Bethe Ansatz 
results obtained in case of infinitely strong Coulomb re¬ 
pulsion [28, 34]. 

Finally, we analyzed the magnetic field dependence of 
the linear conductance and the spectral function, both in 
fields parallel and perpendicular to the CNT axis. Along 
the lines of Ref. [17], we derived the necessary conditions 
to observe the Kondo revival, while extending our analy¬ 
sis to finite valley mixing values (Akk’ 7 ^ 0) and to the 
case of perpendicular magnetic fields. We showed that in 
both parallel and perpendicular fields, the total spectral 
function of the system shows a four peak structure (two 
for the absorption and two for the emission processes). 
We showed that in the case of odd occupancy of a CNT 
shell, a large enough Akk’ can prevent the occurrence of 
a Kondo revival in parallel field; in perpendicular fields 
no revival is expected. Kondo revivals have indeed so far 
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been observed experimentally only for parallel magnetic 
fields [15, 18]. In case of a perpendicular magnetic field, 
the outer peaks can merge with those of the split Kondo 
resonance, thereby leading to a two-peak structure (one 
for absorption and one for emission processes) for large 
magnetic held values compared to the SU(4)-Kondo tem¬ 
perature and B± > Bkk’/2- 
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Appendix A. Orthogonal transformations 


In this Appendix we present the construction of the 
Hamiltonian introduced in Eq. (1) in the Kramers basis 
from an underlying Anderson model that takes into ac¬ 
count the carbon nanotube structure [13, 18]. In the ab¬ 
sence of the valley mixing and spin-orbit interaction, the 
basis set {|.K’ , ,t), , |A,j'), \K, J,}}, indexed by the 

valley and spin quantum numbers, is orthogonal. When 
A ^ 0 this is no longer true, and it is suitable to 
adopt instead the bonding (anti-bonding) representation 
{[a, f), |6,t), |a,4-), 16,4-)}, which can be constructed as 


,a) = |0) = (A.la) 


IM> = dl |0) = 


I K',a) - | K, a) 

V2 


(A.lb) 


In this basis the CNT Hamiltonian consists of several 
terms 


-ffcNT — #CNT + He- e + H B 

= H d + Hkk’ + H S q + H e _ e + H b , (A.2) 

where H d is the SU(4) invariant component and e d the 
orbital energy which can be tuned through the applied 
gate voltage. Explicitly, 


H d = £d 55 • (A.3) 

2,£T=± 

We add to the pure SU(4) term, respectively, the valley 
mixing (Eq. (A. 4a)) and the SOI (Eq. (A.4b)) components: 

H K k- = 55 (A.4a) 

i,<r=± 

H S o = (A.4b) 


Notice that in the bonding/anti-bonding basis the valley 
mixing effect translates in an energy difference between 
the bonding and anti-bonding states [13]. On the other 
hand, the SOI (due to curvature effects and relativistic 
correction to the CNT Hamiltonian) is off-diagonal. 

The fourth term in Eq. (A.2) describes the electron- 
electron interaction with U associated to the charging en¬ 
ergy of the dot 

^ e ~ e = ~2 55 55 4,A<4i,aA,ai- (A.5) 

i=± 

The external magnetic field enters through the last term 
in Eq. (A. 2)). 


Appendix A.l. Spectrum for zero magnetic field 

The single particle Hamiltonian Eq. (A.2) can be di¬ 
agonalized by a unitary transformation, and the new or¬ 
thogonal basis is the Kramers basis {11) , |4), |2) , |3)} used 
throughout of the present work. The unitary operator is 


/ cos ( 6 ) sin (0)0 0 \ 

— sin (0) cos (0) 0 0 

0 0 cos (0) — sin (0) ’ 

y 0 0 sin (6) cos (0) ) 


(A.6) 


where the angle 0 is given by tan(20 ) = Aso/Akk’- In the 
diagonalized form, 77 cnt becomes the dot Hamiltonian 
from Eq. (1). 

We immediately notice that, because the two blocks 
in Eq. (A.6) yield the same eigenvalues, two pairs of de¬ 
generate doublets arise. These are the so-called Kramer 
pairs that, in our notation, are the couples of states (1, 2) 
and (3,4). As such, the states within each Kramers pair 
are related through the time reversal operator T as it 
is sketched in Fig. 2. Additionally, valley reversal, gov¬ 
erned by the anti-unitary operator V, relates the couples 
(1,4) and (2,3) originating from the two sub-blocks. In 
the Kramers basis T and V are given by 

T = k (d\d i + ^4^3 ~ h.c.j , (A.7a) 

V = n (d\d\ — d\d,2 — h.c.^ , (A.7b) 


where k stands for the complex conjugation. Finally, defin¬ 
ing C = V ■ T -1 , it is possible to relate the couples (1,3) 
and (2,4) to each other through the unitary operator 


C = d\d ,3 + d\d 4 + h.c. (A.8) 


Appendix A.2. Spectrum for finite parallel magnetic field 
An external magnetic field parallel to the CNT axis cou¬ 
ples to both the spin degree of freedom and the “orbital” 
one. Thus, in the bonding(anti-bonding) basis, Hb reads 
[18] 

H B =Hl_ s + Hl_ olh 

= 2^ s ^ll 55 a di, crdi,<J + fforb^y 55 ^_ ia di^, 

i,cr—± i, tr=± 

(A.9) 
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where B\\ is the amplitude of the parallel magnetic field. 
Thus, the diagonalized single particle Hamiltonian is 

4 

(B\\) (H & t + #b) U (H||) = Y, AA ( A ' 1Q ) 

4=1 

where 

sin (0 + ) 0 0 \ 

cos (0 + ) 0 0 

0 cos (6~) — sin(0 - ) 

0 sin(0 _ ) cos(6»“) J 

Here we defined OF- in such a way that tan(20 ± ) = (Ago ± 
2g 0 rb-S||)/AKK’- Notice that U [B ||) posses a similar block 
structure as in Eq. (A.6). However, due to 0 + ^ d~, time 
reversal symmetry is broken. 


W(H||) = 


/ cos (0 + ) 
— sin ( 6 + ) 
0 

V o 


Appendix A.3. Spectrum for finite perpendicular magnetic 
field 

An external magnetic field perpendicular to the CNT 
axis couples, differently to the parallel case, only to the 
spin degree of freedom [13, 18]. Its action reads 

Hb = \dsBx Y (A.11) 

i,<7=± 

The transformation that diagonalizes the single particle 
Hamiltonian can be decomposed as a product of two or¬ 
thogonal matrices: 


U(B ± ) = Q 1 (B ± )0 2 (B ± ), 


Or {B±) 


O 2 (B ± ) 


/ cos (0+) 

sin 

(0+) 

— sin 

(0+) 

cos 

(0+) 


0 



0 

V 

0 



0 


u 

0 

1 

°\ 

1 

0 

1 

0 

-1 

71 

1 

0 

-1 

0 


7 

1 

0 

!/ 


0 

0 

cos ( 9 ~) 
— sin ( 6 ~) 


° \ 

0 

sin (9~) ’ 

cos (0~) J 


where 6 ± is such that tan(20 ± ) = Aso/(Akk’ ± g s B±). 
We notice that U (B±) mixes both the spin and the 
bonding(anti-bonding) degrees of freedom, whereas, in the 
previous case, U (-By) was mixing only the latter one. 


Appendix B. Universality and the fixed points 

We devote this Appendix to a brief analysis of the prop¬ 
erties of the SU(2) and SU(4) limiting cases in the N g = 1 
valley. When A = 0 we expect the system to be at the 
SU(4) fixed point. In Fig. B.16 we compare on one side, 
the conductance as computed with the help of Eq. (8) for 
A = 0, with the universal curve for G(T/Tk) as obtained 
from an underlying Kondo model with SU(4) symmetry, 
and we find a perfect agreement. Increasing A, the sys¬ 
tem flows away from the SU(4) and towards the SU(2) 



Figure B.16: Comparison between the temperature dependence of 
the conductance in the limits when A —> 0 and A —> oo (symbols) 
and the SU(2) and SU(4) universal curves for the conductance (solid 
lines). 


fixed point. As the other curve in Fig. B.16 shows, for 
large enough A’s, i.e. A W, the system has already 
reached the SU(2) fixed point. 

Recently, in Ref. [10] a heuristical analytical expression 
was proposed for the universal conductance, of the form 

< bi > 

in order to reproduce the leading-order in the temperature 
expansion predicted by the conformal field theory [35, 36] 
for the SU(4)-Kondo Hamiltonian. In this case the leading 
order is predicted to be cubic, G(T/T K )/G 0 ~ 1 — a ■ 
(T /Tk) 3 + ..., with a a constant of the order ~ 1, even 
though the system has still a Fermi liquid character. This 
sets n = 3 in Eq. (B.l) in contrast to the case n = 2 for 
the SU(2) symmetry. 

In Fig. B.17 we compare the conductance with this 
heuristic curve for different gate voltages. From the fits of 
the DM-NRG data (Fig. B.17), in the range 0 < T/Tk < 1, 
we obtained a good estimate for s = 0.202±0.002 in agree¬ 
ment with Ref. [10]. The heuristic curve reproduces very 
well the DM-NRG results in a wide range of temperatures, 
0 < T/Tk < 10 , and deviations become visible only for 




F g =8.ir — 
K=2.7r — 
Vg= or — 
F g =-2.7r — 

s fit =0.2 


V g =-8.ir - 

T K (V g =8. ir )=0.00 1 6 r 

G 0 (V g =8.ir)=2 \ 

fit - - - 

T K (V g =2.7r)=0.00053r 

G 0 (Vg=2.7r)=2.03 \ 


T K (V r g =0)=0.00043r 

G„(V g =0)= 2 .05 \ 


T K (Vg=-2.7 r )=0.000 4 3 r 

G 0 (V g =-2.7r)=2.07 


T K (V r g =-8. ir )=0.0008 4r 

G 0 (V g =-8.ir)=2.11 



0.001 0.01 0.1 1 10 100 
T/T K (V g ) 


Figure B.17: Universal scaling of the conductance as a function of 
T/Tk for different gate voltages and comparison with against the 
heuristical curve from Eq. (B.l) with n = 3. The other parameters 
are: U = W, T = C//50, T = 0, £d — —U/2 and A = 0. Here Go is 
the unitary conductance in unit e 2 /h. 
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large enough temperatures. 

Trying to fit the data for A W in Fig. B.16 with 
n = 2 we found s = 0.215 ± 0.002, in agreement with the 
predictions from the conformal field theory. 

Appendix C. Numerical renormalization group 
approach 

We solve the Hamiltonian (4) using the numerical renor¬ 
malization group approach. The core of the NRG is a 
logarithmic discretization of the conduction band with a 
parameter A ~ 2, followed by a mapping of the Hamilto¬ 
nian to a semi-infinite chain. In this way the problem can 
be solved perturbatively, as the hopping couplings along 
the chain, P n ~ A _n,/2 decrease exponentially [27]. The 
conduction band Hamiltonian becomes 

4 +oo 

-^chain — + h.c.. (C.l) 

j= 1 n—0 

Here fj n are the fermionic creation operators at the n-th 
site in the j-th channel. The impurity is sitting at site -1 
and is coupled to the site n = 0 

4 

-fftun = V ^ fj,odj + h.c. (C.2) 

i=i 

As the dot Hamiltonian (1) is not modified by this pro¬ 
cedure, the total Hamiltonian consists now of four spin¬ 
less conduction bands coupled to a complex impurity com¬ 
posed of the dot degrees of freedom. 

Along the NRG procedure it is crucial to use the symme¬ 
tries of the system in order to achieve numerically reliable 
results. When A = 0, the model is SU(4) invariant, as the 
total Hamiltonian commutes with the SU(4)-spin operator 

i sw> = [EE W« < c -3) 

n—— 1 

where A = (Ai,..., A 15 ) is a set of matrices defining the 
generators for the SU(4) algebra (generalized Gell Mann 
matrices [37] for example). As discussed in Sec. 2.2, there 
are also two U(l) symmetries corresponding to the conser¬ 
vation of charge in each of the Kramers channels. In the 
NRG language the generators are 

= E (C.4) 

Z n =—1 ' ' 

Since the generators for spin and for the charges com¬ 
mute among themselves, the system has a global U(l) ® 
U(l) ® SU(4) symmetry. The spin orbit and valley mix¬ 
ing perturbations, i.e. A / 0, break the global symmetry 
down to the U(l) ® U(l) (g> SU(2) ® SU(2), generated by 
the charge and the SU(2)-spin operators 

+00 

J* = 2 E E finin' (C.5) 

n=-lj,j'£K 


acting on the two Kramers doublets. Here a = ( cr x , cr y ,cr z ) 
is the vector of Pauli matrices. 
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